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Abstract. We study the interspecies scattering properties of ultracold Li-Cs 
mixtures in their two energetically lowest spin channels in the magnetic field range 
between 800 G and 1000 G. Close to two broad Feshbach resonances we create weakly 
bound LiCs dimers by radio-frequency association and measure the dependence of the 
binding energy on the external magnetic field strength. Based on the binding energies 
and complementary atom loss spectroscopy of three other Li-Cs s-wave Feshbach 
resonances we construct precise molecular singlet and triplet electronic ground state 
potentials using a coupled-channels calculation. We extract the Li-Cs interspecies 
scattering length as a function of the external field and obtain almost a ten-fold 
improvement in the precision of the values for the pole positions and widths of the 
s-wave Li-Cs Feshbach resonances as compared to our previous work [Pires et al.^ Phys. 
Rev. Lett. 112 , 250404 (2014)]. We discuss implications on the Efimov scenario and 
the universal geometric scaling for LiCsCs trimers. 
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1. Introduction 

Universality in few-body systems has been one of the major topics in the ultracold 
quantum gases for the last decade PEI 13]. Its success can be ascribed to the existence 
of magnetic interparticle scattering resonances, called Feshbach resonances (FR), at 
which a two-body bound state crosses the scattering threshold. These resonances are 
routinely employed to tune the interaction strength between the colliding particles P 
and to produce weakly bound dimers by ramping up or down the external magnetic held 
[S]. They can be used to explore intriguing topics in few-body physics, for example, the 
realization of Ehmov’s scenario na Its hallmark is the existence of a geometrical 
series of weakly bound three-body states that exhibits the universal scaling law. The 
energy of the next bound trimer can be found by multiplying the binding energy of the 
previous one with a constant factor. However, the ability to produce and study these 
trimers as well as underlying universal principles, on which the behavior of such exotic 
systems is based, relies on precise knowledge of the properties of the particular FR. 

The central quantity that governs an ultracold collision process, and therefore most 
of the physics at such temperatures, is the two-body s-wave scattering length a. The 
inelastic three-body scattering rate near a FR scales as o^, resulting in magnetic field 
dependent atom losses that can be used to map out how a depends on the external field. 
During the last decade, atom loss spectroscopy in combination with theoretical models 
has become a standard tool in the field of ultracold gases [1| . These methods can give an 
excellent representation of the FR spectrum, however typically not all of the observed 
losses can be unambiguously attributed to an increasing two-body scattering length. 
Especially when a becomes large, not only immediate loss of three atoms from the trap, 
but also other processes, for example, weakly-bound dimer formation and subsequent 
secondary losses may occur. This may lead to shifts and asymmetric broadening of the 
loss signals PEIEIIS] and thus weakening the relation to the functional dependence of 
the scattering length alone. 

More accurate mapping can be obtained by going further than a simple atom loss 
spectrum. The most precise scattering length measurements up to date can be obtained 
by direct radio-frequency (rf) [101 EH (El EHl EH ES] and magnetic field modulation 
13 E3 EH ESI ESI 123 El] spectroscopy of the least-bound molecular states. Since a ER 
intrinsically originates from the coupling of the scattering channel with such a molecular 
state, its energy E in the vicinity of the ER can be connected to the scattering length 
through the relation E ~ a-2 [3122]. By mapping the magnetic field dependence of the 
binding energy of this state, it is possible to study exclusively the two-body problem, 
and the extraction of a is less prone to systematic effects. 

Here we explore the universal behaviour of weakly bound LiCs dimers and Efimov 
trimers close to Li-Cs FRs. We start by investigating the interspecies scattering 
properties of ultracold Li-Cs mixtures of the two energetically lowest LiCs spin channels 
in the magnetic field range between 800 G and 1000 G. We employ rf association and 
atom loss spectroscopy to precisely measure the positions of the LiGs s-wave ERs in this 
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magnetic field range. Depending on the width of the resonance we separate them into 
two groups and use complementary approaches to determine their properties. For the 
broad resonances close to 843 G and 889 G we measure magnetic field dependent binding 
energies of weakly bound dimers through rf association. For the narrow resonances we 
employ atom loss spectroscopy. Due to their small width the resonance position can still 
be detected with high accuracy. We use these measurements as an input for a coupled- 
channels (cc) calculation that allows us to construct accurate Li-Cs molecular potentials, 
from which scattering lengths, resonance positions and widths are determined. The 
obtained parameters agree well with the previous observations [231 121] and recent 
extensive theoretical studies [23], however they represent almost an order of magnitude 
improvement in precision and accuracy. Finally, we discuss the implications of these 
results to the recent observation of LiCsCs Efimov resonances [2Sl EZj • With the help 
of the new mapping of the scattering length in the vicinity of the 843 G LiGs FR we 
obtain refined Efimov scaling factors of 5.5(2) and 5.0(1.5) for the first and second 
Efimov period, respectively, where the first period slightly deviates from the universal 
Efimov scenario of 4.9, as predicted for the Li-Cs system with mass ratio of 22 [2H|- 

2. Radio-frequency association of LiCs Feshbach dimers 

The sample preparation scheme for the rf association measurements is similar to the one 
presented previously in Refs. [231 [26] . In brief, we prepare the Li-Cs mixture in a crossed 
optical dipole trap with standard laser-cooling techniques. Using degenerate Raman 
sideband cooling [23] most of the Cs atoms are optically pumped and spin-polarized in 
the energetically lowest spin state |F = 3,mF = 3). After the initial cooling steps the 
Li atoms populate both energetically lowest spin states, namely |F=l/2,mi7’ = l/2) 
and |F = 1/2, = —1/2). The last forced evaporation ramp is performed at 920 G. 

At the end of the ramp one of the two Li spin components is selected by shining in a 
short, resonant light pulse that expels the other one from the trap. Finally, about 3 • 10^ 
(4-10^) atoms remain in the respective Cs (Li) spin channels with a temperature around 
400 nK for each species. We measure the trapping frequencies /a,, fy, fz of 11 Hz, 114 Hz, 
123 Hz (33 Hz, 275 Hz, 308 Hz) for Cs (Li) atoms, where the external magnetic field is 
parallel to the z axis. The total uncertainty of the applied magnetic field amounts to 
16 mG (one standard deviation), resulting from long-term magnetic field drifts, residual 
field curvature along the long axis of the cigar-shaped trap and calibration uncertainties. 

To associate the molecules we start with a mixture prepared in the non-resonant 
scattering channel at a variable magnetic field close to the FR in the resonant state (see 
Fig. [^. We drive the system with a rectangular rf pulse with the frequency Erf/h that 
is close to the resonance frequency Eo/h between the two energetically lowest Li spin 
states. In order to determine the molecular binding energy Ef,, we scan the frequency of 
the applied rf field and observe the number of Li atoms that are left in the non-resonant 
state after the rf pulse. 

A typical loss spectrum is depicted in Fig. Detuned from the free-free transition, 
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Figure 1. Radio-frequency association of Li-Cs molecules. The mixture is initially 
prepared in the non-resonant scattering channel, here Li|l/2,—1/2) 0 Cs|3,3), at a 
magnetic field close to the broad 843 G s-wave Feshbach resonance in the resonant 
scattering channel Li |l/2, l/2)0Cs |3, 3), which couples to the weakly bound molecular 
state under study. Depending on the frequency of the rf driving field either free- 
free (with the energy Eq) or free-bound (with the energy Eq 0 Ej,) transition can 
be studied. An analogous scenario is implemented close to the second broad Li-Cs 
Feshbach resonance in the Li |l/2, —1/2) 0 Cs |3,3) scattering channel close to 889 G. 
In this case the mixture is initially prepared in the Li |l/2,1/2) 0 Cs |3,3) channel. 


which corresponds to a flip of Li nuclear spin, we observe an additional loss feature that 
originates from the association of LiCs Feshbach dimers (free-bound transition). We 
also identify a similar loss signal at comparable values of detuning and amplitude in the 
remaining number of Cs atoms. To limit saturation effects, we experimentally optimize 
the power and length of the association pulse such that at most 30% of atoms are lost 
at the end of the rf pulse. The optimized pulse length ranges from 0.5 s close to the FR 
and 7 s away from it. 

3. Theoretical analysis 

3.1. LiCs dimer association spectrum 

We model the observed loss spectrum with the help of rate equations and the dimer 
binding energy dependent two-body association rate [301 [ini. The long association 
pulse lengths and low molecule yield, which is below our detection limit, indicates that 
the dimer association rate is much smaller than their loss rate. Assuming a quasi¬ 
stationary state, in which each produced molecule immediately gets lost through atom- 
dimer collisions, the time dependent Li atom loss at a given magnetic held is governed 
by and thus can be described through the following equation 

Nu = 


( 1 ) 
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Figure 2. Remaining number of Li atoms after an rf association (free-bound) pulse of 
LiCs molecules at a magnetic field of 842.04 G and relative Li-Cs temperature of 400 
nK. Each data point is an average of three measurements and the error bars represent 
the standard error. The binding energy is determined from the fit of Eq. to the 
data (solid line) and yields = 69.7(1.6) kHz and 7 = 5(1) kHz for a pulse length 
t = 3 s. The vertical dashed line corresponds to the fitted binding energy. 


where is the initial number of Li atoms in the non-resonant state, rics denotes the 
density of the Cs gas cloud, t is the length of the applied rf pulse and contains 
the functional form (see below) of the molecular association rate. Here we neglect 
single-body losses and assume constant rics^ which is justified by the fact that we are 
working in a low saturation limit of the rf transition. By solving the full system of rate 
equations we estimate that this approximation, using the simple Eq. [T} may introduce 
a minor error on the fitted atom loss amplitude, which does not exceed 10 %. Since the 
number of produced LiCs molecules at any given point through the experimental cycle 
is insignificant, we do not include the loss terms associated with the molecule-molecule 
recombination. 

The two-body association rate is determined by the energy-dependent wave- 
function overlap of the scattering atom pair with the final molecular state [30]. For a 
thermal ensemble it can be expressed as 


i^rf) ~ ^ j h{Sr)F(Sr^ Eq + + Sr)dSr^ 

0 


( 2 ) 


where h{sr) oc e is the number density of colliding atom pairs with relative energy 

Sr and temperature T, and 


E{Sr-) E/fe) C)C 


V^rEbEl 

V^y {Sr + EbyiSr + E') 


( 3 ) 


is the energy normalized Franck-Condon density between the scattering wave function 
of a free Li-Cs atom pair and a bound Feshbach dimer with binding energy Eb [HOI [ID]. 
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E'l^ is defined through the Li-Cs reduced mass fi and the non-resonant channel scattering 
length a' as E[ = The convolution of the spectroscopic line shape with the 

Lorentzian profile L^{Erf, Eq + Ef, + Sr) of width 7 accounts for the strong collisional 
broadening, yielding an estimated lifetime of LiCs molecules in the mixture around 
30 fis. The prefactor C contains all the numerical factors resulting from the integration 
of rate equations, and experimental parameters that affect the molecule production rate, 
but which are approximately constant for a given magnetic field, as well as species- 
dependent atom-dimer inelastic collision rate^ It also accounts for uncertainties in 
the determination of the absolute gas densities, which, under realistic experimental 
conditions, can vary up to a factor of two due to systematic errors in measurements of 
the trap frequencies, temperature and the exact number of atoms. 

The binding energy of the Feshbach dimers at a given magnetic field is extracted 
by fitting Eq. to the loss spectrum of Li atoms, as displayed in Fig. We use 
Eb, 7 , and C as free fitting parameters and set a' = —28.5 Uq [2Sl EH ES]- 
Small variations in a' that are of the order of a few of percent have affect the fitted 
binding energies on a permille level. The temperature of each species is determined 
in an independent measurement with identical trapping parameters and is kept fixed 
during the fit. To exclude systematic effects associated with the precise determination 
of relative temperature we verify that by increasing it by a factor of two the value of 
the fitted binding energy does not change by more than 1 kHz. By performing the 
measurements and the fitting procedure for different external magnetic fields, we record 
the binding energy dependence, which is displayed in Fig. for the two broadest FRs 
in the Li-Cs mixture. 

The extracted binding energy can be affected by several other systematic effects. 
One of them is the mean-field shift, which starts to dominate in the regime where the 
scattering length is comparable to the interparticle spacing, i.e. na^ ~ 1. For our 
experimental densities of n ~ 10 ^^ cm“^ such shifts would become relevant at magnetic 
field regions with binding energies on the order of Eb ~ 1 kHz, which is sufficiently far 
away from the region where the experiments were performed. Additionally, by changing 
the background Cs atom density we checked that the observed molecular association 
line shifts stay within the statistical uncertainties of the fit, and therefore we do not 
include the mean-field shift in the analysis. 

Another source of systematic resonance shifts is the confining optical dipole 
potential. The detuning of the dipole trap laser beam is large in comparison to the 
hyperfine splitting of the involved spin states, hence its created fight shift is equal for 
both of them and can be neglected. However, the confining potential can contribute to 
the scattering state energy shift in two other ways. The first one is the confinement 
induced shift of the relative ground state energy for two colliding free atoms. Its 
magnitude can be calculated for two interacting particles in a cigar shaped harmonic 
trap bh, and for our dipole trap geometry with aspect ratio rj ^ 9 it yields 325 Hz. The 

§ The prefactor C depends on ncs and For atom losses, which do not exceed ~ 30% of the initial 
number of atoms, it does not change by more than 5 %. 
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Figure 3. Binding energies of LiCs Feshbach molecules and atom losses. The left 
and right panel correspond to the magnetic field regions near the 843 G and 889 G 
Feshbach resonance, respectively. The blue crosses display the dimer binding energy 
that is extracted from a fit of Eq. 0to the rf association spectrum at the given 
magnetic field. The error bars represent one standard deviation of the total error, 
which results from statistical and systematic uncertainties. The blue and red solid 
lines show the calculated molecular state energies from the coupled-channels model 
and the universal binding energy Ei, = with the resonance parameters 

from Table respectively. The red shaded region corresponds to the uncertainty 
of the FR parameters. The green squares show the remaining Gs atom number for a 
corresponding atom loss measurement. Here the error corresponds to one standard 
error of the mean. The systematic magnetic field uncertainty for the atom loss 
measurements in this figure is 30 mG. The vertical dashed line displays the resonance 
pole position, and the gray shaded region corresponds to the uncertainty. The arrow 
in the left panel shows the position of the second excited LiGsGs Efimov resonance 
with scattering length a_ ’. 


second complication is the fact that, in general, the problem of two different atoms in 
a harmonic trap with unlike trapping frequencies does not separate into center-of-mass 
and relative coordinates. The magnitude of the shift of the associated lowest energy 
state in a Li-Cs mixture can be estimated for our trapping geometry and mass-ratio, 
and is on the order of 50 Hz j32]. Since the order of magnitude of these corrections is 
much smaller than the measured binding energies we neglect these effects in the model 
that we use to fit the data, however, we include them in the total systematic error 
budget. 

To obtain a complete set of Li-Cs s-wave FR properties we re-measure the positions 
of the narrow s-wave resonances in the two energetically lowest scattering channels up to 
1000 G. For these measurements we reach roughly an order of magnitude lower relative 
kinetic energies than in the previous works [231 Eg. Their experimental positions R® 
are determined with a Gaussian fit to each of the loss features in the magnetic field 














range where the line shape is approximately symmetric. The improved magnetic field 
stability and lower temperatures allows us to determine these resonance positions with a 
roughly five-fold better precision than in our previous measurement [23]. These results 
are summarized in Table [H 

3.2. Coupled-channels caleulation 

To obtain an accurate mapping between the scattering length and the magnetic field that 
is independent of the employed analytical fitting model [3S] , we analyze the data with the 
full cc calculation for the Li(2s)+Cs(6s) asymptote as in the previous work [2HIE5|- In 
short, the determination of the final resonance positions relies on the creation of accurate 
LiCs molecular potential curves for the electronic singlet and triplet a^E'*' ground 

states. The potentials are constructed in such a way that they simultaneously reproduce 
the binding energies of the Feshbach molecules, the refined s-wave FR positions 5® from 
atom loss experiments, and 6498 rovibrational transitions from laser-induced Fourier- 
transform spectroscopy [31|. We deduce the theoretical resonance positions from 
the maxima of calculated two-body collision rates at the experimental kinetic energy. 
For the binding energies below the 843 G resonance we exclude the two data points 
with the smallest binding energies from the fit. Their rf association spectra, due to 
experimental limitations, already overlap with the Li free-free transition spectra, which 
hinders a reliable extraction of the free-bound spectra (similar to the one in Fig. for 
these respective magnetic field values. 

The results of the modeling are listed in Table [I] as deviation 5 = B'^ — B* from 
the measured positions for the experimentally employed relative collision temperature T 
and drawn as solid lines in Fig. These results provide almost an order of magnitude 
improvement over the previous determination of the FR positions through the trap-loss 
measurements [23l |2l| and the rf spectroscopy [26], and they are consistent with the 
recent theoretical analysis [2H], if the differences in determining the resonance positions 
and experimental accuracy are taken into account. 

Finally, we characterize the resonance profiles by calculating the scattering length 
dependence on the magnetic field at a kinetic energy of 1 nK and fitting this dependence 
with the conventional functional form 

“ ^ “ B-Bfr “ ■ ■ ■) 

with as many terms as there are resonances in the given channel. The resonance 
position Bfr, its width A, and background scattering length Obg are used as free fitting 
parameters, and they are given in Table [T] By including all observed resonances in a 
single fit of Eq. one removes a possible slope of the effective background scattering 
length resulting from neighboring resonances. We note that one could also fit the 
calculated profiles by a product of resonance functions for each resonance instead of the 
sum. This will result mainly in different values of A, but as long as one is using only the 
derived functional values, the interpretation is consistent. The fitted values reproduce 
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Table 1. Positions of the LiCs s-wave Feshbach resonances. Unless specifically 
noted, the experimentally obtained resonance positions are extracted by fitting a 
Gaussian profile to the loss spectra for the relative collision temperature T, at which 
the measurements were made. The numbers in the brackets represent the total error 
that includes uncertainty of the magnetic field, and statistical and systematic errors 
of determining the position of the resonance. The results of the coupled-channels 
calculation are given as deviations 6 = B^ — B^ from the observations and show 
excellent agreement with the data. Bfr, A, and a^jg gives the fitted resonance pole 
position, width, and background scattering length, respectively, for the calculation 
with kinetic energy of 1 nK. 


Entrance channel 

(G) 

5(G) 

T (nK) 

Bfr (G) 

A(G) 

dbg (do) 

Li|l/2,-l/2) 

816.128(20) 

-0.005 

300 

816.113 

-0.37 

-29.6 

0 Gs |3, H“3) 

888.595(16)* 

0.002 

100 

888.578 

-57.45 

-29.6 


943.020(50) 

-0.033 

400 

943.033 

-4.22 

-29.6 

Li|l/2,+l/2) 

842.845(16)* 

-0.000 

100 

842.829 

-58.21 

-29.4 

0 Gs |3, 03) 

892.655(30)t 

0.005 

100 

892.629 

-4.55 

-29.4 


* Extrapolated from rf association. The temperature shown is only used for the 
calculation of the scattering resonance and selected sufficiently low to reduce its 
influence to less than 5 mG. The error reflects the uncertainty of the field calibration. 
^ This measurement was performed in a double-wavelength optical dipole trap with 
species selective optical potentials. Details will be given elsewhere. 


the calculated s-wave scattering length to better than 2% in the entire magnetic field 
range between 500 G and 1000 G, which we also use for the re-evaluation of the Efimov 
resonance positions (see the next section). There is a slight difference between the two 
theoretically obtained resonance pole positions and Bfr. We strongly suspect that 
it originates from the different types of numerical calculations that were used to extract 
these parameters, however further investigation is necessary to find the exact reason 
behind this difference. Therefore we estimate the error for the resonance pole positions 
from the systematic error from the magnetic field calibration and difference between the 
theoretical values, which yields ± 23 mG. 

We use the resonance parameters that were obtained from the cc calculation to 
plot the simple single-channel formula Eh = fi^/(2/ia^), which relates the universal 
dimer binding energy to the scattering length (see Fig. |^. The Li-Gs characteristic van 
der Waals energy scale jl] is 157 MHz, thus the influence of the short range effects on 
the measured binding energies is minimal. This is reflected in the nearly ideal 1/a^ 
scaling of the measured binding energies in this magnetic field range. Since the Li- 
Gs background scattering length ahg ~ — 29.5ao is small and negative, we expect only 
very minor influence of the virtual state in this regime. This contributes to the simple 
situation where the two LiGs binding energies are well described with the universal 
relation and can be treated independently from other neighboring resonances in the 
same scattering channels. This is in contrast to more complicated situations, like the 
one in Gs atoms where FRs overlap [IHl IS5] • 
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Table 2. Positions and scaling factors of the LiCsCs Efimov resonances near the 
843 G Feshbach resonance. The error of the magnetic field, at which the resonances 
are found, represent the combined statistical and systematical error arising from the 
determination of the Efimov resonance positions and magnetic field uncertainty. For 
and the scaling factors the uncertainty of the a{B) mapping through the Feshbach 
resonance parameters is added to the error. 


Efimov state i 

Magnetic field (G) 

a^- (ao) 

(i) , (i-1) 

a_la_ ^ 

0 (ground) 

848.90(7)t 

-311(4) 

- 

1 

843.85(3)t 

-1710(70) 

5.5(2) 

2 

843.03(6)t 

-8540(2650) 

5.0(1.5) 

0 

848.55(12)* 

-329(7) 

- 

1 

843.82(5)* 

-1760(100) 

5.3(3) 

2 

842.97(4)* 

-12200G140) 

6.9(2.3) 


t Ref. [26] 
* Ref. HZ] 


Our determined position of the 843 G FR pole clearly deviates from the previously 
observed atom loss maximum [2S1 EZ], as can be seen in Fig. It also deviates from 
the result BpR = 842.75(3) G obtained by Tung et al. [27] . where exclusively atom loss 
measurements are used to infer the resonance pole position. This illustrates that the use 
of atom loss alone is questionable for a reliable determination of the FR pole position, 
especially if the resonance width is much larger than experimental uncertainties, as it 
is in the present case. The dehnition of the resonance pole position is based on pure 
two-body scattering whereas a number of different loss mechanisms may contribute to 
the total loss effect, the most prominent being the three-body collisions. The situation 
in the vicinity of the resonance’s pole is complicated furthermore by the fact that not all 
of inelastic three-body collisions result in an immediate loss of the atoms from the trap. 
Contribution from other recombination processes, for example, weakly-bound dimer 
formation and subsequent atom-dimer recombination, should be considered, which may 
lead to increased loss away from the pole of the FR [H| . In this case the maximum of total 
atom losses can be shifted with respect to the maximum in the corresponding three- 
body loss rate. The specific loss channels and the exact pathway of this decay in Li-Cs 
system still remain an open question, requiring a selective product state determination, 
which is not available at the present stage of our experiment. However, we expect that 
the shifts between the determined scattering pole positions and experimentally observed 
loss maxima can be explained or influenced by similar mechanisms as those discussed 
for other systems of ultracold gas mixtures pmiHii. 

4. Scaling of LiCsCs Efimov resonances 

With the improved determination of the 843 G FR pole position we re-evaluate the 
scattering lengths and scaling factors of the previously observed Efimov resonances 
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[261 EZ], In Table we summarize the reported magnetic field values and corresponding 
scattering lengths at which the z—th three-body Efimov state merge with the 
scattering threshold. The extracted scaling factors of 5.5(2) and 5.3(3) for the first 
Efimov period are close to the expected value of 4.9 for a zero-temperature gas in 
the scaling limit (|a| ^ a) pfl EH], however they slightly deviate from the universal 
prediction. This is not surprising, since the assumption of the scaling limit is not strictly 
justified for the ground state Efimov resonance [2S|- Several theoretical studies have 
demonstrated that the ground Efimov state can be subject to large modifications due 
to short range physics [HHl EZl EHl ESI SO] and even three-body forces m- Eurt her more, 
the recent study of the first excited-state resonance in Cs [121SS] and the new analysis 
of ®Li data mi do not only give a hint to deviations from the universal scaling, but 
also to shifts of the ground-state resonances due to finite range effects. The exact origin 
of the above mentioned deviation in the Li-Cs system still remains an open question. 
The present analysis also does not consider Effmov resonance shifts arising from ffnite 
temperature. This is in contrast to the recent studies in the homonuclear gases of Cs and 
Li atoms mum SHI , which have shown that such effects need to be taken into account 
in order to accurately determinate the Effmov resonance positions. At the same time, 
the scaling factor of the second Efimov period lies well within the universal predictions, 
and experimental demonstration of deviations from the universal law, if any, will require 
access to lower temperatures and improved control over the magnetic field systematics. 

5. Conclusion 

In conclusion, we have precisely determined the resonance positions of all the s-wave ERs 
in the two energetically lowest Li-Cs scattering channels. The present work represents 
almost an order of magnitude improvement in accuracy and precision over the previous 
determination. It was achieved by performing rf association of colliding Li-Cs atom 
pairs into weakly bound dimers close to the broad ERs around 843 G and 889 G and 
complementary atom loss measurements of further three narrow s-wave resonances. 
Based on the measured magnetic field dependent binding energies and atom losses 
precise singlet and triplet molecular potential curves for the LiCs electronic ground 
state were constructed with the help of a cc calculation. The obtained potentials were 
used to map the Li-Cs scattering length on the external magnetic field. The precise 
resonance parameters will be pivotal for future experiments in the Li-Cs systems, where 
they can serve as a starting point for the preparation and investigation of ultracold 
polar molecules [ffilllTlIlH], the creation of strongly interacting mixtures and polarons 
[12l[HniEIllS2], and enable access to the Efimov physics in the universal regime [211121 ESI- 
Additionally, the accurate knowledge of the Li-Cs ER parameters has allowed us 
to re-evaluate the positions of previously observed Efimov resonances. We have found 
a slight deviation from the predicted universal scaling factor for the first period of 
LiCsCs Efimov trimers, while the scaling factor of the second period is consistent 
with the universal law. This result is intriguing, since it approaches a regime where 
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the applicability of few-body theories to a realistic mass-imbalanced system can be 
tested quantitatively. To fully enter it, however, lower temperatures will be necessary. 
Presently observed Efimov features may be significantly infiuenced by the short range 
interactions and the unitary limit, which both obscure clear observation of the second 
excited LiCsCs Efimov resonance and might have strong effect already on the first 
excited one. This may lead to shifts of the resonance positions and require more 
sophisticated models for accurate description of the three- and more-body physics near 
FRs. 
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